Jet energy loss in the quark-gluon plasma by stream instabilities 
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We study the evolution of the plasma instabilities induced by two jets of particles propagating in 
opposite directions and crossing a thermally equilibrated non-Abelian plasma. In order to simplify 
the analysis we assume that the two jets of partons can be described with uniform distribution 
functions in coordinate space and by Gaussian distribution functions in momentum space. We find 
that while crossing the quark-gluon plasma, the jets of particles excite unstable chromomagnetic 
and chromoelectric modes. These fields interact with the particles (or hard modes) of the plasma 
inducing the production of currents; thus, the energy lost by the jets is absorbed by both the gauge 
fields and the hard modes of the plasma. We compare the outcome of the numerical simulations 
with the analytical calculation performed assuming that the jets of particles can be described by 
a tsunami-like distribution function. We find qualitative and semi-quantitative agreement between 
the results obtained with the two methods. 

PACS numbers: 12.38.-t,12.38.Mh 



I. INTRODUCTION 

The properties of the matter produced in heavy-ion collisions can be studied in several different ways [l[ , by observing 
electromagnetic signals, by studying the dissociation of heavy quarkonia states {QQ) or by studying the propagation 
of high pt partons. The high px partons generated by hard scatterings in the initial stage of the heavy- ion collision 
behave as hard probes of the produced medium [2J. When a jet of partons travels across the medium it loses energy 
and degrades. The study of the mechanism of energy loss gives interesting information about the properties of the 
traversed medium (see [3[ for reviews). The process of absorption of the energy of the jet by the quark-gluon plasma 
(QGP) has been named jet quenching, because in the direction of propagation of the jet one observes a decrease of 
high energy hadrons, whilst the number of soft hadrons increases. 

Several models have been proposed to describe the jet quenching mechanism, see Q for a review and detailed 
discussions. In previous papers 0, 0] we have presented a novel mechanism for describing how the jet loses energy 
and momentum while traveling in a thermally equilibrated quark-gluon plasma. Since the jet of particles is not in 
thermal equilibrium with the QGP it perturbs and destabilizes the system inducing the generation of gauge fields. 
Some of these gauge modes are unstable and grow exponentially fast in time, absorbing the kinetic energy of the jet. 
Plasma instabilities in the QGP have been the subject of extensive work in the recent years, see Q for a review. This 
is so because it was suggested 0, HI that the presence of plasma instabilities could help to explain fast equilibration 
of the QGP. 

In Ref. we have studied the propagation of an energetic jet of particles in an equilibrated quark-gluon plasma by 
using both transport theory and fluid equations [j| . In both cases we have considered many simplifying assumptions. 
We have neglected hard modes interactions, assuming that the interaction between the jet and the plasma is only 
mediated by mean gauge fields. Moreover, the plasma has been considered to be in thermal equilibrium, whereas the 
jet is assumed to be initially electrically neutral and colorless and described by a uniform distribution function in 
space and a tsunami-like distribution in momentum space. 

In the present paper we consider the same simplifying assumptions reported above in a setting where two jets 
of particles propagate in opposite directions in a thermally equilibrated non-Abelian plasma. We first analytically 
derive the growth rate of the unstable gauge field modes and then compare the results with the outcome of numerical 
simulations. In the numerical simulations we describe the particles of the jet employing a uniform distribution function 
in space and a Gaussian distribution function in momentum space. By varying the width of the Gaussian distribution 
function we can approximate the tsunami-like distribution function and then compare the results with the outcomes 
of the analytical study. We find qualitative and semi-quantitative agreement between the results of the two methods. 
Moreover, the numerical study gives some hints regarding the dynamical processes responsible of the energy loss of 
the jets. We find that the hard modes of the jet induce the instability of the gauge fields of the plasma, meaning 
that the energy associated with the chromomagnetic and chromoelectric fields increases exponentially fast in time. In 
turn, the increase in time of these modes induces currents of the particles of the plasma. Therefore, the energy lost 
by the jets is absorbed by both the gauge fields and the particles of the plasma. This is a remarkable result, since in 
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our model there arc no direct interactions (collisions) among the hard modes of the jets and the hard modes of the 
plasma. The only interaction that is present is the one mediated by the gauge fields. 

It is worth recalling that the interaction of relativistic streams of particles with (abelian) plasmas has been studied in 
different fields of physics, ranging from inertial confinement fusion, astrophysics and cosmology. In these experiments 
one observes that when the particles of the stream are charged, plasma instabilities develop, leading to an initial stage 
of fast growth of the electromagnetic fields. These instabilities have been studied using a variety of methods, from 
kinetic theory to hydrodynamics [Tol. [Tl|; experimental evidence of the relativistic filamentation instability has been 
reported in Ref. [12j |. 

This paper is organized as follows. In Section [TT] we describe the model and review the basic set of equations. In 
Section [IIII wc derive the analytical solutions for the case of two jets described by a tsunami- like distribution function 
that travel across the equilibrated quark-gluon plasma. Moreover, we derive some general results regarding the energy 
density of hard-particles. In Section HVl wc present the results of the numerical simulations. We draw our conclusions 
in Section fVl 



II. FLUCTUATIONS IN THE QGP TRAVERSED BY A NEUTRAL JET OF PARTICLES 

For the sake of simplicity, we consider a plasma consisting only of gluons (adding quarks and anti-quarks is straight- 
forward, as shown in that is initially in a colorless state described by an isotropic distribution function. (The 
treatment of anisotropic plasmas can be done as in [Hi ], where one assumes that the momentum of particles is 
symmetric only for rotation around an axis.) 

For a thermally equilibrated and isotropic plasma we use the Bosc-Einstcin distribution: 

/pi = IbeM = eP J T _ 1 > (1) 

and a uniform distribution function in coordinate space. 

Without loss of generality the particles that constitute the jets are taken to be masslcss quarks that are initially in 
a colorless state described by a uniform distribution function in coordinate space and a distribution function /j c t(p) 
in momentum space. Several distribution functions have been proposed [Til ] typically expressed in the power-law form 

where po is introduced as an infrared cutoff and n is some positive number. In order to facilitate the analysis of the 
dynamical properties of the system in Ref. Q a tsunami-like distribution function was used 

f iet (p)=nS^(p-A) , (3) 

that describes a system of particles of constant density n, all with the same momentum p = A. In the analytical 
study of the gauge field instabilities we shall describe the jets by means of this tsunami-like distribution function. 

In the numerical study we use a smooth function of momentum, which in a certain limit becomes a tsunami-like 
distribution function. In this way we can compare the results of the numerical simulations with the analytical results. 
More specifically, we use the Gaussian distribution 

/£ t (p)=iV(a)e-^- A > 2 , (4) 

where N(a) is a normalization factor and a measures the smearing in momentum space. In the limit a — > oo, one has 
that jjo t (p) — > /jet(p) and all the particles of the jet have the same momentum A. We fix the normalization factor 
by requiring that the particle density is the same for both the tsunami-like and the Gaussian distribution functions, 
namely 

which is fulfilled for 

N(a) = n(47ra) 3/2 . (6) 

We shall numerically study the system evolution for various values of a and n, in order to determine under which 
conditions the plasma instabilities are stronger. It can be ex pec ted from general arguments that plasma instabilities 
are present for any non- vanishing value of a, sec e.g. Ref. [15| . however it is important to determine whether the 
growth rate is sufficiently large so that these instabilities can play a significant role in the evolution of the system. 
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A. Deviation from equilibrium 



The two jets of particles that propagate in the non-Abelian plasma excite unstable chromomagnetic and chro- 
moclcctric fields destabilizing the system. In this Section we derive the equations governing the small deviations 
from equilibrium of the distribution functions of gluons and hard modes. The approach we present below describes 
successfully the physics of the soft scales in the plasma, or, the Hard Thermal Loop (HTL) physics as discussed in 
Ref. [IB, E3|- We extend that formalism to add the contribution of the two jets of particles. 

Small deviations from equilibrium of the distribution functions of the gluons of the plasma can be described by an 
(N% — 1) x (JVg — 1) matrix in color space SG(p, x), where N c is the number of colors. For the particles of the jet the 
deviation from equilibrium can be described by a N c x N c matrix in color space SWj a t(p,x). 

The linearized transport equations for the gluons and the particles of the jet can be written as 

P^V^Gip, X) = -gp^ v ( X )^L , p^SWjetip, X) - -gp»F^(x)^ , (7) 

where the covariant derivatives T> ^ and act as 

V ll =d ll - ig[A^(x), ... ] , D^=d^- ig[A^(x), ... ] , 

with A M = A^(x)r a and A^ = A%(x)T a , and r a and T a are the SU(N C ) generators in the fundamental and adjoint 
representations, respectively. The strength tensor in the fundamental representation is F^ u , while T^v denotes the 
field strength tensor in the adjoint representation. 

For massless particles, it is possible to do the following rc-definitions 

5G(p,x) = -g^W£{v,x) , 6W ist ( P> x) = -g^W» t (v,x), (8) 

then considering that the distribution function of the plasma is isotropic one obtains that the transport equations for 
the gluons of the plasma in the previous equations simplify to 

p^SG(p,x) = -g^ { X) ^M } SG( P ,X) = -g^M W O i(VjX) ; (9) 

apo apo 

and therefore W£(v,x) = 5 vQ W^{v,x). 

Regarding the jet, it is characterized by an anisotropic distribution in momentum space, therefore one has that both 
W® ct (v, x) and W? et (v, x) arc non vanishing. Upon substituting the expressions of Eqs. §8§ in Eqs. J7J one obtains the 
following equations 

v*VpW£ = v^ v , v^D^ = , (10) 

where w M = p^/p with p = |p|. 

Since the fluctuations of the distribution function of the plasma and of the jet arc not in general color singlets, it 
follows that one has to construct the color currents associated with these fluctuations 



8%{x) = -g Jp» r a Tr[T a SG(p,x)] , Sj^(x) = 8W iet (p,x) , (11) 

where the momentum measure for massless particles is defined as 

The color currents generated by the fluctuations of the distribution functions of the plasma and of the jet enter as 
sources in the Yang-Mills equation 

D,F^ (x) = 6% (x) = 6f p (x) + 6j? et (x) . (13) 

Equations (fT3"|) and Eqs. (TIT)]) form a set of equations that has to be solved self-consistcntly. The gauge fields which 
are solutions of the Yang-Mills equation enter into the transport equations of every particle species and, in turn, affect 
the evolution of the distribution functions. 
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III. ANALYTICAL SOLUTION FOR THE TSUNAMI-LIKE DISTRIBUTION FUNCTION 



Analytical solutions for the color polarization tensor and for the corresponding dispersion laws of the collective 
modes can be obtained when the jet is described by the tsunami-like distribution function in Eq. ([3]). The analysis of 
the system constituted by a single jet crossing an equilibrated plasma has been done in Ref. [f|. Here we extend that 
calculation to the case of two anti-parallel jets crossing the plasma. 

We briefly recall some results of the HTL approach and derive the dispersion laws that describe the evolution of 
gauge collective modes [l8| . The dielectric tensor is defined as 



Tl ij 



(14) 



where in order to simplify the notation, we have dropped the color indices. For the QGP in the HTL approximation 
one has that 



e «( w> k) = (a* 



k l W 



k l W 

e T (w,k) + -p-£ L (w,k) , 



where 



£L(w,k) 
e T (w,k) 



2k 2 



uj cu — k 
2k ° S u) + k 



lu k \ cu — k 

log 

2k 2uj uj + k 



(15) 

(16) 
(17) 



are the longitudinal and transverse (with respect to k) components of the dielectric tensor of the plasma composed 
by gluons and massless quark of Np flavors, and ui 2 = ^g 2 T 2 (N c + Np/2) is the plasma frequency squared. 

Now we consider the collective modes which are present in a system described by the tsunami-like distribution 
function ([3]). In this case one deduces that for the two jets propagating in opposite directions one has that 



n& x (fc) = 



fcV + fcV {J 2 - k 2 )v z v j 
~ (w - k • v) 2 



n« a (fc) = -wf M [ 5 lJ - 



"jet 



uj — k • v 

fcV + kHi 1 



(cj 2 



k- v 



k- vf 



(18) 
(19) 



where the plasma frequency of the jet is given by uj 2 ct = ng 2 / (2A) and the two jets propagate with the same velocity 
|v| in opposite directions. 

In the analysis of the unstable collective modes we are interested in very short time scales when the Vlasov 
approximation can be employed. The effect of the beam of particles is to induce a color current, which provides a 
contribution to the polarization tensor. Therefore, apart from the HTL polarization tensor, we have to include the 
polarization tensor associated to the jets. For very short time scales, the polarization tensor of the system is given by 
the sum of the polarization tensors of its components, meaning that 



= (fc) + nj^fc) + uf: t2 (k) ; 



thus, the dielectric tensor of the system is given by 

e ij {u;,k) =6 ij 

and the dispersion laws of the collective modes can be determined solving the equation 

det Ik 2 ^' - k i y - oj 2 e ij (k)\ = . 
The solutions of this equation depend on |k|, |v|, cos 9 = k ■ v, on 

2 2 9 

u: = u„ + ui ' 



P "jet 



(20) 



(21) 



(22) 



(23) 



and on 



~jet 



(24) 
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We consider two different cases corresponding to k || v and k _L v and determine the behavior of the unstable gauge 
modes of the system solving Eq. . The results for the imaginary part of the dispersion law as a function of k for 
the case k |j v arc given in the left panel of Fig[TJ We consider various values of v and find that the largest instability 
is obtained for v ~ 0.9. The corresponding results for the case kl v are reported in the right panel of Fig[T] Notice 
that in both cases the range of validity of our results is k < Lu t , but we have also reported the results for larger 
momenta. In the case k _L v it happens that for large values of k the growth rate saturates; however, the saturation 
is achieved in a range of values of the momenta where the approximations that we are employing are not reliable. 
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k/GX 



FIG. 1: (color online) Imaginary part of the dispersion law of the unstable modes of the system composed by two jets of 
particles propagating in opposite directions in a non-Abelian plasma. The various curves are obtained solving Eq. (|22p for 
various values of the velocity and taking b = 0.01, see Eq. (|24fl . Left panel: longitudinal component, k || v. Right panel: 
transverse component, k _L v. The full (black) line corresponds to v = 0.5, the dotted (red) line corresponds to v = 0.7, the 
dashed (green) line corresponds to v = 0.9, and the dot-dashed (blue) line corresponds to v = 0.999. These solutions for the 
unstable modes have been obtained assuming that the jets of particles have a uniform distribution in space and the tsunami-like 
distribution function in momentum space reported in Eq. ©. 



In Table [I] we report the largest value of the growth rate of the longitudinal and transverse components for various 
values of b. Both the longitudinal and transverse growth rates increase with increasing values of b. The transverse 
modes always have a growth rate that is larger than the one of the longitudinal modes. 



b 


F||M 


i±M 


0.01 


0.02 


0.12 


0.1 


0.05 


0.4 


0.2 


0.06 


0.6 


0.5 


0.08 


0.9 



TABLE I: Approximate largest value of the imaginary part of the dispersion law for the parallel and transverse cases, for various 
values of b. For a given value of 6 the largest value of Fy corresponds to v ~ 0.8, whereas for T± it corresponds to i) = 1. 



A. Energy density of the hard modes 

The contribution of the quasiparticles (hard modes) to the energy-momentum tensor of the plasma can be deter- 
mined using the following equation 

Tr(Sj flt Fn=9^ + d fl Qf: t , (25) 

where Sj^t is the variation of the total current defined in Eq. (|13[) . and and ®j et are the energy-momentum 
tensors of the plasma and of the jet, respectively. Then, the procedure proposed by Blaizot and Iancu in Ref. [20l | can 
be used to single out the energy-momentum tensor for a plasma close to equilibrium, and express the result in terms 
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of W®i- The same strategy can be generalized for treating anisotropic, non-equilibrium plasmas as shown in Ref. [l9l ]. 
For the evaluation of the energy density of the plasma one can use the equation 

Trfe t i^) = 4/^3 {^ W > ' D M + f» ' D )K] > (26) 
and since we are assuming that the plasma is isotropic, we have that 

Trfe^°) = 4 / Ifa {^gf W > ' + • ■ (27) 

Integrating by parts the first term in the equation above one can identify the energy density associated to the plasma 
obtaining, in agreement with the results of Blaizot-Iancu, that 



g 2 f d 3 p df p i 
(2tt) 3 d Po 

In this linear analysis, it is clear that 8™ is always positive. However, in Eq. (|28[) the fluctuations of the gauge 
fields have been neglected. When these fluctuations are sufficiently large they can give the dominant contribution to 
Opj* and eventually can make this quantity negative. This is indeed what we have observed in the early stage of the 
numerical simulations that will be presented in the next Section. However, for sufficiently large times the contribution 
of fluctuations is not dominant and 0°° is positive in agreement with Eq. (|28|) . Regarding the energy density of the 
jet it is not possible to obtain a close expression analogous to Eq. (f2"5| . Therefore, in order to establish whether 
the energy density associated to the jet is positive or negative one has to rely on numerical simulations. On general 
grounds one expects that if the quasiparticlcs of the jets are sufficiently energetic they will lose energy while crossing 
the plasma. The point is to observe whether the rate at which the jet will transfer energy to the plasma will be linear 
or exponential with time. In the second case it means that the jet has been able to trigger an instability of the system. 

IV. NUMERICAL SIMULATIONS 

In this Section we present the numerical simulations of the real-time evolution of the system composed by the 
two jets of particles propagating in an equilibrated and isotropic plasma. For the sake of simplicity we consider the 
non-Abelian SU(2) gauge group and restrict the gauge field dynamics to 2 + I dimensions. 

A detailed description of the numerical algorithm that we employ can be found in Ref. [l9| . In order to implement 
the numerical simulation we introduce the quantities 

, 9 2 r 00 dp 2 df 

av = "Wo ^ p dp-> (29) 
g 2 f 00 d P 2 df 



^ = -W„ SF^' (30) 

for both the plasma and the jet particles. Since we are considering an isotropic plasma, one obtains that 6^ = 0, 
while a^ 1 = m 2 D , where mu is the Debye mass. 

The two jets propagating in opposite directions are characterized by the Gaussian distribution function of Eq. ((U) 
with Ai = (0, 0, A) for the first jet, while A2 = (0, 0, —A) for the other. For the first jet we find that 



i) 



^2 e -a(p 2 +A 2 -2p, 2 A) ^ (32) 



that can be explicitly evaluated obtaining 



AT p / 2 \ 

^ = ^ ^erf(L) V 1 + ^ + yzV * eV * (3/2 + + ° rf {Vz)) ) ' (33) 

= -^ ^-^ (y. + ^(1/2 + vl)(l + erf(y,))) , (34) 
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FIG. 2: (color online) Plot of ccjj jt / w jet (left panel) and of fo^V^jet (right panel) defined in Eqs. (|33[1 as a function of v z for four 
different values of L. With increasing values of L both functions become more and more peaked at v~ — 1. 



where erf(x) is the error function, where 

L = V^A , (35) 

is a dimcnsionlcss variable, y z = v z L and oj? ct — 4^ is the plasma frequency of the jet. (This definition of the plasma 
frequency is the same we have used in the analytical study.) Notice that a v and b v depend on w? et and on functions 
of L of the form L m cxp(— L 2 ), where m — 1, 2. The color current associated with the jet can now be expressed as 

Kt M = / £ v " (<*Kt («. x ) + K et W e t («,*)), (36) 

where the angular integral can be numerically evaluated by using the discretization of the unit sphere with 

Zi = -l + {2i-l)/N z , i = l...N z ; tff^tocj/Nt, j = l...N$. (37) 

The resulting "disco balls" are such that they are covered with AT = x N z tiles of equal area. In order to test 
the robustness of the algorithm we have performed the integrals increasing the values of and N z , until numerical 
stability is achieved. 

A plot of the functions and 6^* is presented in Fig.[5]for four values of L. Both a^* and are peaked at v z = 1 
and the value at the peak increases with increasing L. Moreover, the region where a^* and fe^* arc non vanishing 
becomes narrower with increasing values of L and for L —* oo they are non zero only for v z = 1 where they diverge. 
For this reason one cannot numerically evaluate the integral in Eq. (|36[) for arbitrary large values of L. Therefore when 
doing the simulation one cannot consider large vales of L; we shall take L < 6. 

Notice that 6^* < for any value of L and v z . For L = one has that = w? et , while fa^' vanishes. The reason is 
that for L = there is no flux of particles, the distribution function of the jet becomes symmetric. The corresponding 
expressions for the second jet can be obtained replacing A — > —A in the equations above. 

An expression similar to Eq. (|36[) holds for the current of the plasma, however the evaluation of the integral is 
simpler because we have assumed that the plasma is isotropic. Notice that the anisotropies in the system are due to 
the two jets, but the system is still symmetric for rotations around the direction of propagation of the jets, which 
move in opposite directions. This means that one can study the gauge field dynamics in 2 + 1 dimensions, with the 
two space dimension orthogonal to the direction of propagation of the jet. This allows to reduce the simulation time. 



1. Results of the numerical simulation 



The simulation depends on several parameters. In Fig. [3] we present a typical plot of the time evolution of the 
absolute value of the energy of the chromomagnctic fields, of the chromoelcctric fields, of the harddoop modes of the 
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FIG. 3: (color online) Time evolution of the absolute value of the energy of longitudinal and transverse chromomagnetic fields, 
B z and B t , of the longitudinal and transverse chromoelectric fields, E z and E t , and of the hard- loop modes of the plasma, 
HL-pl, and hard-loop modes of the jet, HL-jet. In this simulation we have used N z = 30 and N$ = 30, see Eq. (|37[1 . and the 
parameters b = 0.2 and L = 4, see Eqs. p4[) and (|35[) . The full black line corresponds to AE, that is the energy violation 
measured during the evolution of the system. The results of the numerical simulations become unreliable when AE becomes 
of the same order of the total energy of the system. 



plasma and of the hard-loop modes of the jet. The black line corresponds to the violation of the energy measured 
during the evolution of the system. 

The most unstable gauge field mode corresponds to the transverse chromomagnetic component. This is in agreement 
with the analytical results we obtained in Section IIIII where it was found that the strongest instability appears in 
gauge field modes with momentum orthogonal to the jet velocity. 

However, we can now look more in detail to what happens in the system. The hard modes of the jet lose energy 
that is absorbed by the gauge fields of the plasma, that exponentially grow. Notice that the growth of the gauge 
fields of the plasma drives the growth of the hard modes of the plasma. This effect is due to the interaction between 
the gauge fields and the hard modes of the plasma. In order to understand this point one has to consider that in 
our simplified model there are no collisions between particles, therefore hard particles of the plasma do not directly 
interact with the hard modes of the jets. The only interaction between the two hard modes is mediated by the gauge 
fields of the plasma. In other words, what the hard modes of the plasma feel is an increase of the magnetic field and 
respond by generating a current and absorbing energy from the gauge fields. 

We can schematically say that the sequence of events is the following: The hard modes of the jet induce the 
instability of the gauge fields of the plasma, meaning that a chromomagnetic field is induced that increases in time. 
Then, the increase in time of the chromomagnetic field induces a current of hard modes of the particles of the plasma. 

This sequence of events can be deduced/confirmed looking at Fig. [31 where one sees that the energy associated to 
the hard modes of the plasma starts to grow only when the energy of the gauge fields is sufficiently large. 

A peculiar effect can be observed in the early stage of the evolution. Contrary to what happens at late times, 
for very short times the hard modes of the plasma lose energy, whereas the hard mode of the jet gain energy. This 
behavior is due to the initial fluctuations of the hard modes of the plasma and of the jet (that induce fluctuations of 
the gauge fields as well). However, this effect can be modified changing the initial conditions of the system. As an 
example if we put fluctuations to zero we find that the energy associated with the hard modes of the plasma is always 
positive in agreement with Eq. (|28[) . 

Regarding the quantitative analysis, the growth rate of the gauge fields depends on b and L. In Table HT1 we report 
the values obtained for various values of these quantities. 

The results of the numerical simulation are in qualitative and semi-quantitative agreement with the analytical 
results. In both cases the growth rate of the most unstable modes correspond to the transverse ones and the growth 
rate increases with increasing values of b. Moreover, the values of the growth rate in the two approaches are quite 
similar as can be seen comparing Table U with Table [Til This is a remarkable result, because in the numerical 
simulation one has an integrated growth rate, meaning that modes with all momentum compatible with the lattice 
cutoff contribute to the growth rate. On the other hand, the growth rate of the analytical case has been determined 
considering the largest value of the curve in Fig. [1] 
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b 


L = 2 


L = 4 


L = 6 


0.1 


0.09 


0.21 


0.24 


0.2 


0.12 


0.29 


0.31 


0.5 


0.15 


0.32 


0.35 



TABLE II: Values of the growth rate, T/ut, of the unstable gauge fields for various values of b and L. 



V. CONCLUSIONS 



We have studied a system composed by an equilibrated and isotropic quark-gluon plasma traversed by two jets of 
particles propagating in opposite directions. We have analyzed the instabilities of this system in two different ways. 
First, we have carried out an analytical study using a very simplified description of the distribution function of the 
jets, the so-called tsunami-like distribution, while the plasma has been described by a Bose distribution function. 
Then, we have performed a numerical study, assuming that the distribution functions of the jets are two Gaussian 
distribution in momentum. In both cases we have assumed that the plasma and the jets have uniform distribution 
in coordinate space. In order to simplify the numerical analysis we have restricted the gauge field dynamics to 2 + 1 
dimensions and we have considered the SU(2) gauge symmetry. 

Comparing the results of the two approaches we have found qualitative and semi-quantitative agreement for the 
behavior of the unstable gauge fields. In both cases the most unstable modes correspond to the transverse chromo- 
electric and chromomagnetic fields. Moreover, in both approaches we find that the growth rate is about an order of 
magnitude smaller than the total plasma frequency uj t , defined Eq. (|23p . and increases with increasing values of &, see 
Eq. ([24]) . The plasma instabilities fully develop on time scales of the order of t ~ 10 fu>t and we can estimate the upper 
bound of this time scale evaluating the plasma frequency in a weak coupling scenario. Considering a temperature 
T ~ 300 MeV we find that t ~ 1 fm/c. 

From the numerical simulations one realizes that the energy lost by the hard modes of the jets is absorbed by both 
the gauge fields and the hard modes of the plasma. It happens that the hard modes of the jet induce the instability 
of the gauge fields of the plasma and therefore the energy associated with the chromomagnetic and chromoelectric 
fields of the plasma increases exponentially fast in time. Then, the increasing chromomagnetic and chromoelectric 
fields induce currents of the particles of the plasma and transfer energy to these hard modes. 

The analysis we have performed is simplistic in many ways. In particular, we have assumed that both the plasma 
and the jets have a uniform distribution in space. It should be possible to improve the treatment presented here in 
future studies. One should localize the jets in space-time, with the probable associated effect of having a localized 
energy deposition. So far, we have only considered the very short time evolution of the system. At sufficiently long 
times, radiative and collisional processes should be included, with the addition of collision terms in the transport 
equation. In a weak coupling scenario, the collisional times arc suppressed by a factor 1/g 2 with respect to the time of 
development of the instabilities. Collisions will probably stop the growth of the gauge fields, but they will represent 
a later mechanism of jet energy loss. It might be interesting to see if one could incorporate in a numerical code all 
the above mentioned effects, along the line presented in Ref. [2l[ 

Even with the simplifications employed in this paper we believe that the effect we discuss, which is based on 
collective plasma phenomena, might be relevant for the description of the energy loss mechanism of relativistic jets 
traversing a quark-gluon plasma. 
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